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Overheating is a failure mode that significantly affects the reliability of 
electronic devices. All electronic devices, including a 3-phase inverter 
driving a traction motor, produce heat dissipation. Heat dissipation needs to 
be controlled with cooling to prevent overheating. Overheating can be 


avoided by increasing cooling or reducing heat dissipation. Heat dissipation 
in the 3-phase inverter is caused by the internal resistance of the metal— 
oxide-semiconductor field-effect transistor (MOSFET), switching loss, and 
Keywords: other factors. Cooling for the 3-phase inverter can use water coolant or air 
coolant. The cooling system is based on the amount of heat dissipation 
produced. Cooling of a 3-phase inverter can use air coolant with the addition 
of an extended surface area in the heat sink. The heat sink uses aluminum 
Extended surface material, often called pin fin. There are kinds of aluminum available in the 
Heat dissipation market. We calculated heat generation based on the MOSFET's internal 
Heat sinks resistance, switching loss, and other factors. We validated the simulation 
results experimentally using a thermal camera. Thus, we could find an 
optimal number, dimensions, and aluminum type of fin for the cooling 
system in the 3-phase inverter. 
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1. INTRODUCTION 

Electronic component reliability is influenced by a cooling system [1]. Each electronic component 
will produce heat dissipation, including the 3-phase inverter for a motor controller. Heat dissipation needs to 
be controlled with a cooling system to prevent overheating. It can also be done by reducing the amount of 
heat dissipation produced, such as using the right switching strategy [2]-[4], using field-oriented control 
(FOC) [5], [6], using intelligent control [7]-[11], and the others. The most dominant heat source in the motor 
controller is losses in the 3-phase inverter. Part of the 3-phase inverter that will produce heat is the metal— 
oxide-semiconductor field-effect transistor (MOSFET). 

One of the MOSFET failure modes is overheating. The heat generated by MOSFETs must be 
transferred quickly to the cooling media with intelligent cooling using a fan [12], thermoelectric cooling [13], 
others so that the temperature of the MOSFET does not rise above its ability limit. Heat transfer can occur in 
several ways, including the following: conduction, convection, and radiation [14]-[16]. This paper will 
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discuss the MOSFET cooling method using an extended surface (pin fin). The extended surface area of heat 
transfer improves the thermal efficiency [17], [18]. 

Computational fluid dynamic (CFD) simulation can be used to obtain a cooling method that suits your 
needs [19], [20]. This method can help to analyze the cooling method without spending time doing 
experiments [21]. The numerical method can show the temperature distribution and how much heat must be 
cooled [22]. If the final design has been obtained, this design can be used as a reference in making a prototype. 
This work presents optimized heat transfer from MOSFET in a 3-phase inverter by using extended surfaces with 
some pin fins and dimensions. The calculation of heat loss in the 3-phase inverter is used as the input in the 
simulation. The simulation result is validated with the experimental results. The validation is carried out based 
on the temperature distribution and the temperature's value at the same point. If the error tolerance that occurs 
does not exceed 596, we consider that an accurate simulation is obtained. Therefore, it can be used to optimize 
the number of fins and the dimensions of the fins. Our target is that the 3-phase inverter can operate normally at 
a power of 30 kW. 


2. METHOD 

First, we performed a heat loss calculation in 3-phase inverter for parameter input in the simulation. 
Next, the simulation results were validated with the experimental results, where the temperature distribution 
was obtained using a thermal camera. After the simulation results are validated, it is possible to optimize the 
heat sink according to the heat loss. As a result, the optimal heat sink configuration is obtained with the 
controller not overheating. 

Heat losses can be seen in the power flow diagram in Figure 1. Based on Figure 1, the 3-phase 
inverter losses consist of conduction loss pulse-code modulation PCM, switching loss part submission 
warrant PSW, dead time loss PD, and gate charge loss PG, as given in (1) to (4) [23], [24]. 


Pew = ze fo?" poy(t)dt = =f)" Roson: Edt = Roson: rms (1) 
Pow = <Viw-lo. (ty + ty). few (2) 
Py = Vp. lo. (tor + tof). few (3) 
Pg = (Qg-u F 0): Vos- fsw (4) 


The symbol Rpson is the MOSFET drain-source internal resistance (Q), Tsw is switching time (s), ipis the 
current flowing in drain (A), and /p,4; is the current flowing from source to drain (A). Vin and I, represent the 
battery voltage (V) and the output current (A). t, is the rise time of the MOSFET (s), tris the fall time of the 
MOSFET (s), and fsw is the switching frequency (Hz). Vp, tpr, and tp; denote the low-side MOSFET — diode 
forward-voltage (Volt), dead-time at rising (s), and dead-time at falling (s). Vgs is the gate drive voltage of the 
MOSFET (V), Q; uis the high side MOSFET gate electric charge (C), and Q,. ; is the low side MOSFET 
gate electric charge (C). This controller uses the type of MOSFET IRFB4115 with the specifications shown 
in Figure 2. 


Base Metal 
| Rpsíon) max :11 mQ 


Maximum Temperature : 175 °C 
Turn-On Delay (taon) : 18 ns 
Rise Time (t;) : 73 ns 
Turn-off Delay (tao) | : 41 ns 


Deadtime Loss Plastic Body 
Conduction Loss 
Switching Loss 


Power Input Gate-Drive Loss Power Output 
Pa = I, Š Va Pout = 3. Inmsa - Vrmsa -cos PHI 
Figure 1. Power flow 3-phase inverter Figure 2. MOSFET IRFB4115 specification [25] 
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This controller has the following target specifications; i) maximum voltage operation 150 VDC, ii) 
maximum power 30 KW, and iii) using an air-cooling system. We first made a simplified 3D model of the 3- 
phase inverter for numerical analysis using CFD software, as shown in Figure 3. The heat source on the 3- 
phase inverter is the MOSFET component. So that heat transfer can be analyzed with thermal resistance as 
shown in Figure 4. 


Thermal Pad 


MOSFET x 48 
Base Plat 


Fin Heatsink 


Figure 3. Simplification of 3D model of the 3-phase inverter 
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Figure 4. Thermal resistance model controller 


The heat source from the MOSFET moves either up to the controller cover or down to the heat sink. 
The heat transfer downward will be more dominant than upward. The controller cover traps the air contained 
inside the controller container, so it does not allow air circulation. The heat transmitted downward propagates 
through the thermal pad. 

A thermal pad is a component that serves as an insulator between the base metal of the MOSFET 
and the heat sink. The use of the thermal pad must have good thermal conductivity as well as a good 
insulator. The motor controller failure can also be caused by improper thermal pad selection. If the heat 
generated by the MOSFET is not transferred to the heat sink quickly, the MOSFET will experience 
overheating, resulting in failure. Similarly, if the thermal pad cannot isolate the metal base of the MOSFET, a 
short circuit will occur, and the controller will fail. Generally, the thermal pad uses silicone material. The 
heat transfer then continues to the heat sink fin by conduction. The equation of heat transfer by conduction is 
given in [14], [16]: 


Qconduction = k. A. AT (5) 


The type of heatsink material affects the quality of heat transfer: the more significant the heat sink material's 
thermal conductivity (k), the more remarkable its heat propagation ability. Our heat sink uses 1060 series 
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aluminum material, often found in the local market. Convection sends the heat transfer into the air around the 
heat sink fin. The equation of heat transfer by convection is given by (6) to (8) [14], [16]: 


Qconvection = h. A. AT (6) 

Nu = = = f (Re,, Pr) (7) 
.V.L 

Re, == (8) 


The symbols h and A represent the convection coefficient (W/m?K) and the heat transfer area (m°). L 
is the length of material (meters), ky is the thermal conductivity of the fluid (W/mK), and Nu is the Nusselt 
number. Re and Pr denote Reynold's number and Prandtl’s number. p is the density of air (kg/m?), V is the 
velocity of air (m/s), and u is the friction coefficient. 

Convection heat transfer can be optimized by increasing the convection coefficient and cooling area. 
Increasing the convection coefficient can be done by increasing the speed of air flowing at the fin heat sink, 
for example, by adding a fan, directing the air with a guide fan so that air can be directed precisely at the 
heatsink. The cooling area can be enlarged by using variations in dimension and the number of fins on the 
heatsink. The quality of the simulation results is determined by the dimensions and meshing used. Our 
meshing uses a curvature-based parameter mesh with total nodes of 615,925, using Jacobian points by 4 
points. Table 1 lists material types used in the motor controller cooling simulation. Figure 5 shows the 
meshing model. 


Table 1. Material types for cooling simulation 


Component Materials Thermal conductivity (W/mK) 
Fin heat sink Aluminum 1060 200 
Base plat Aluminum 1060 200 
Thermal pad Silicone 32 
MOSFET plastic body Plastic Silicone 124 
MOSFET base metal Aluminum 5052 137 


Figure 5. Meshing model 


Connection contact set using bounded type. The contact that needs to be defined is between 
MOSFET with the thermal pad and the thermal pad with the heat sink. This step ensures that the connections 
between different materials are perfectly connected. We conducted numerical simulation using CFD with the 
following conditions: 

— Surrounding temperature: 28 °C 
— Controlled variable: Convection coefficient (h) 
— The source of the heat loss: the MOSFET 

First, we calculated heat loss when the motor controller is under a 1.5 kW load as show in Table 2. 
Afterward, the simulation can be run, and temperature distribution can be obtained. Simulation results will be 
validated using experimental results. Experiments with a maximum load of 30 kW continuously (30 minutes) 
cannot be done because the dynamometer will overheat. So, the experiment was carried out with a continues 
load of 1.5 kW. Furthermore, the results of these experiments can be used to validate the simulation results. 
The validated simulation results can then be used to test the maximum loading, so that the temperature 
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distribution is obtained when the maximum loading to obtain the temperature distribution under the 
maximum loading. This simulation can also be used to optimize the fin on the heat sink so that cooling is 
obtained as needed with minimum fin dimensions. 

Figure 6 shows the experimental setup with (a) schema of experiment and (b) installed equipment in 
this paper. It is equipped with an eddy current dynamometer paired with a 30 KW motor. The temperature 
sensor uses a 10K thermistor recorded in a datalogger (USB to CAN Datalogger). An oscilloscope with four 
channels is used to measure voltage and current, both input and output. A thermal camera is used to obtain 
the temperature distribution. 


Table 2. Calculated heat loss in the motor controller under 1.5 kW load 


Type Heat Loss Parameter Value 
Conduction Loss Roson 11 mQ; Ip= 8.15 A; Vin = 80.34 V; 0.617 W 
Switching Loss t,—73 ns; t¢=39 ns;fsy-20kHz; 0.733 W 
Deadtime Loss fat = 18 ns; taom =41ns; Vgs = 23V; Qg =120 nC. 0.772 W 
Gate Charge Loss 0.055 W 

Total Heat Loss 2.177 W 
Coupling 
Phase 
bu S dicii BLDC Motor | =" 
E 
E 
Current Phase E 
E BLDC Motor Controller Thermal Camera 
o: ill wtape Phase E 
scilloscope [ee ee ( 
‘Current Battery USB to CAN 
Datalogger 
Computer Battery 
Pack 
Computer 
(a) 


Figure 6. Experiment setup (a) schema and (b) equipment 


3. RESULTS AND DISCUSSION 
3.1. Simulation validation through experiment 

The load was kept constant at 1.5 kW with a motor speed of 950 rpm and variations in wind speed. 
We changed wind speed by controlling a fan power to obtain speed 1, speed 2, and speed 3. CAN to USB 
datalogger was used to obtain data, among others: temperature of heat sink, phase current, battery current, 
battery voltage, and the motor speed. One-time data collection took 35 minutes or until the heat sink 
temperature got steady. The result of the experiment is shown in Figure 7. 
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Figure 7. Motor speed and heat sink temperature 
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The heat sink temperatures in Figure 7 were measured using a thermistor sensor. Yellow, red, blue, 
and green line curves are the heat sink temperatures with variations in wind speed 1, speed 2, speed 3, and 
natural convection. Each line's steady-state temperature is yellow at 33.37 °C, red at 33.96 °C, blue at 
34.92 °C, and green at 39.72 °C. The greater wind speed yields a lower temperature of the heat sink. This is 
in accordance with (6), (7), and (8). The trend line that occurs tends to be steady at a specific time. This 
indicates that the cooling is sufficient so that overheating will not occur. The temperature distribution can be 
captured in this steady state with a thermal camera. 

Figures 8 and 9 show the temperature distribution results obtained from the simulation and 
experiment. Figure 8 shows the temperature distribution using h=40 W/m?K, both experimentally and in 
simulation results. Figure 8(a) shows the simulation results where the thermistor sensor laying point shows a 
temperature of 33.1 °C. While Figure 8(b) is a photo taken by the thermal camera where the location of the 
thermistor shows the value of 33.3 °C. These results indicate that the simulation results have a 0.6% error. 
The temperature distribution in general also has similarities. Under continuous loading of 1.5 kW, the highest 
temperature occurs in the DC-DC converter, yielding 37 °C. 


Temp (Celcius) 


33.1 Celcius $2 


(a) (b) 


Figure 8. Temperature distribution with h=40 W/m?K (a) simulation result and (b) experiment result 


Figure 9 shows the temperature distribution using h=33 W/m?K for both experimental and 
simulation results. Figure 9 (a) shows the simulation results where the thermistor sensor laying point shows a 
temperature of 34.3 °C. In comparison, Figure 9 (b) is a photo taken by the thermal camera where the 
location of the thermistor shows the value of 34.9 °C. The result indicates that the simulation results have a 
1.75% error. Based on the results shown in Figures 8 and 9, the simulation results can be accepted with an 
error tolerance of + 1.75%. Furthermore, the simulation method can be used to evaluate the maximum 
loading of the 3-phase inverter, i.e., with a load of 30 kW. Table 3 summarizes the data obtained from the 
simulation and experimental results. 


- Temp (Celcius) 
343 Celcius 36 
| 37.1 


(a) 


Figure 9. Temperature distribution with h=33 W/m?K (a) simulation result and (b) experiment result 


(b) 
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Table 3. Comparison of simulation with experiment results 


Coefficient of Convection Experiment Simulation Percentage of Error 
h= 33 W/ m'K 34.9 °C 34.3 °C 1.75 % 
h= 40 W/ m'K 33.1 °C 33.3 °C 0.6 96 


3.2. Optimization of the heat sink under maximum load 

We calculated heat loss produced by each MOSFET when the 3-phase inverter operates under the 
maximum load 30 kW. Table 4 shows the parameter values and the heat loss values. Since we have 48 
MOSFETs, the total heat loss is 1,048.27 Watts. We used this heat loss value as the input to our CFC model 
and conducted a numerical simulation. Figure 10 shows the temperature distribution of the heat sink under a 
maximum load of 30 kW. The maximum temperature is 98.9 ?C, so the controller does not overheat. The 
temperature where the MOSFET is situated shows a value of 38.5 °C. Since the temperature is significantly 
lower than the critical temperature of the MOSFET, the heat sink fin is not required in this area. 


Table 4. Heat loss of each MOSFET in the motor controller under 30 kW Load 


Type Heat Loss Parameter Value 

Conduction Loss Rosíon= 11 mQ; Ip = 65.16 A; Vin = 80.34 V; 14.491 W 
Switching Loss t,=73 ns; tp=39 ns; fgy=20kHz; 3.551 W 
Deadtime Loss taon) = 18 ns; taom =41ns; Ves = 23V; Qg =120 nC. 3.742 W 
Gate Charge Loss 0.055 W 
Total Heat Loss 21.839 W 
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Figure 10. Simulation result under maximum load 30 kW 


Furthermore, we optimized the number and length of the heat sink fins with maximum thermal load 
and h=60 W/m?K. The original length of ach heat sink fin is 20 mm. Figure 11 shows the simulation results. 
Figure 11(a) shows the temperature distribution using 15 pieces of heat sink fins. The highest temperature of 
the MOSFET is 103 °C. The temperature is lower than the critical operation temperature of the MOSFET, 
which is 125 °C. Figure 11(b) shows the temperature distribution using 12 pieces of heat sink fins. The 
highest temperature of the MOSFET is 107 °C. The result indicates that reducing of fin number to 12 pieces 
yielding maximum temperature that is lower than the critical temperature of the MOSFET. Further 
optimization is carried out by reducing the fin length to 12 mm and yielding the temperature distribution 
shown in Figure 11(c). The maximum temperature of the MOSFET in this condition is 117 °C, which is 
lower than the critical temperature of MOSFET. However, there is an increase in temperature compared to 
the result of 12 pieces 20 mm fin. Therefore, the controller under a maximum load of 30 kW is optimally 
cooled using the configuration of 12 pieces of fins with 12 mm length, following Figure 11(c). 
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Figure 11. Optimized heat sink in the number and length of the fins: (a) 15 pieces of 20 mm long, 
(b) 12 pieces of 20 mm long, and (c) 12 pieces of 12 mm long 
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4. CONCLUSION 

Thermal analysis of the controller cooling system was carried out using modeling and simulation. 
The simulation was validated experimentally with an error of + 1.75%. Furthermore, the simulation method 
can be used to optimize the cooling design for a maximum loading of 30 kW. The results indicate that the 
controller does not experience overheating, as demonstrated by the MOSFET temperature. Further 
optimization, through simulation was carried out in order to reduce the size and the number of the fin. The 
simulation result of the optimized heat sink suggests that the maximum temperature of the MOSFET is lower 
than its critical temperature. Hence, the controller did not experience overheating. 
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